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Abstract. We address the question of quantifying entanglement in pure graph states. 
Evaluation of multipartite entanglement measures is extremely hard for most pure 
quantum states. In this paper we demonstrate how solving one problem in graph 
theory, namely the identification of maximum independent set, allows us to evaluate 
three multipartite entanglement measures for pure graph states. We construct the 
minimal linear decomposition into product states for a large group of pure graph 
states, allowing us to evaluate the Schmidt measure. Furthermore we show that 
computation of distance-like measures such as relative entropy of entanglement and 
geometric measure becomes tractable for these states by explicit construction of closest 
separable and closest product states respectively. We show how these separable states 
can be described using stabiliser formalism as well as PEPs-like construction. Finally 
we discuss the way in which introducing noise to the system can optimally destroy 
entanglement. 
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1. Introduction 

Ever since the realisation of entanglement's importance in quantum information 
processing a large effort has been devoted to classifying states according to their 
entanglement properties pQ. This has proven to be a daunting task since unlike in the 
bipartite case we encounter a much richer structure when characterising multipartite 
quantum states [2j [3]. Even when one concentrates on a coarser picture based on 
separability properties the question remains formidably hard to settle [3]. 

Entanglement quantification in multipartite quantum states is one of the 
fundamental problems in quantum information theory. In this work we concentrate 
on three measures of genuine multipartite entanglement, namely the relative entropy of 
entanglement [5], geometric measure [6j[7J[8], and the Schmidt measure [9]. As all three 
measures are defined as minimisations of distances in Hilbert space or over all linear 
decompositions into product states they are extremely hard to compute analytically. In 
this paper we restrict ourselves to pure graph states. 

Examples of states for which any of these measures can be computed are sparse and 
usually contain some form of symmetry or admit an efficient description that facilitates 
the evaluation. The class of symmetric states is one such example for which the relative 
entropy of entanglement and the geometric measure can be computed [TU1 [LL|- The 
form of the closest separable state for pure and mixed cluster states, a particular regular 
instance of a graph state, has been investigated in [T2J [13]. More general treatments 
of finding the closest separable state were also attempted by inverting the problem 
and asking what the closest entangled state is given a separable state on the boundary 
between entangled and separable states [HI EE] . However these approaches are limited 
to the scenario of two qubits. 

Computation of relative entropy of entanglement and the geometric measure for 
some classes of graph states has been considered in [12]. The Schmidt measure has 
been found for most of 7-qubit pure graph states equivalence classes generated by local 
Clifford transformations in [16]. This approach was extended to 8 qubits in [17J. A 
common feature among these treatments is that the amount of entanglement was found 
indirectly by computing the lower and upper bounds for the respective measures. In 
many cases the bounds coincide and therefore the exact value for entanglement can be 
obtained. 

We extend and modify these techniques to show how the entanglement measures 
can be computed directly. This may seem as a daunting task at first but we demonstrate 
how the optimisation problem of evaluating the entanglement measures can be mapped 
to a well-known problem in graph theory of finding the maximum independent set, 
or equivalently the minimum vertex cover. Doing this allows us to approach the 
optimisation problem from a more graphical and intuitive perspective. Furthermore 
this strategy reveals a previously unrecognised connection between the closest separable 
state and the minimum linear decomposition of the graph state into product states. 

Structure of the paper is the following. In Section 12.11 we review the three 
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multipartite entanglement measures evaluated in this paper and known relationship 
between them. For the reader's convenience we also give a brief discussion of graph 
theory and some quantities that are central to our argument in Section 12.21 Section 12.31 
offers a quick overview of graph states and stabiliser formalism. Finally we conclude 
this review section by a discussion of lower and upper bounds on the three entanglement 
measures in Section \2M 

In Section [3] we present our main results and discuss in detail how the entanglement 
measures can be evaluated directly. As we are able to evaluate the three entanglement 
measures only in the case when the lower and upper bounds coincide we discuss the 
necessary and sufficient conditions when this is the case in Section 13.11 Section 13.21 
presents our main results and describes how to evaluate the three entanglement measures 
directly. The main logic of our argument is to identify a suitable stabiliser subgroup 
of the original graph state stabiliser. This subgroup stabilises a subspace spanned by 
product states that can be used to construct the closest separable state and hence 
evaluate the relative entropy of entanglement. Furthermore this subspace contains the 
original graph state as the equal superposition of the spanning product states allowing 
us to find the minimal linear decomposition of the graph state vector into product 
states and hence to compute the Schmidt measure. Knowing the minimum linear 
decomposition we can easily identify the closest product states that maximise the overlap 
with the pure graph state making it possible to evaluate the geometric measure. 

Even though the stabiliser formalism offers the most comprehensive picture when 
constructing the closest separable state it is limited in the sense that it only applies 
to pure graph states. Therefore in Section H] we explore two other descriptions of the 
closest separable state. The first description is based on a construction similar to that 
of projected entangled pairs (PEPs) [18]. Usually PEPs are used to describe pure 
entangled states whereas we show how the same argument can be used to describe 
mixed separable states. The main motivation behind this approach is to generalise our 
previous techniques to states that cannot be described within the stabiliser formalism 
such as weighted graph states [19]. Our second description of the closest separable states 
considers how one can optimally destroy entanglement by introducing noise to the graph 
state. 

Finally in Section |5] we summarise our work and present possible directions of future 
research and some open questions. 

2. Preliminaries 

In this section we review the three entanglement measures considered in this paper and 
some of their relevant properties. As the derivation of majority of our result relies on 
concepts from graph theory we also briefly review the relevant quantities. Finally we 
summarise the main properties of graph states including how they can be described 
using the stabiliser formalism and how the bounds of entanglement measures can be 
calculated. 
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2.1. Entanglement measures 

Relative entropy of entanglement (REE) of an entangled state p is defined as follows [5] 
E R (p) := mm S{p\\u), 

where S := Tr[plog p — p\ogu\ is the quantum relative entropy and the minimisation is 
taken over all separable states u. Since we investigate only pure states in this paper the 
relative entropy takes a simplified form of S = — Tr[plogo;]. We refer to the separable 
state that achieves the minimum of the relative entropy as the closest separable state 
(CSS). Even though the relative entropy is not a true metric we can interpret E R {p) 
as the shortest distance between the entangled state p and the set of separable states 
SEP. 

Geometric measure [6j [7J E| for a pure state can be defined as 
E a (\1>)):= mm - log |(</#>| 2 , 

\4>)€PROD 

where the minimisation is taken over all product states Note that in the case of 
geometric measure the minimisation is taken over pure states and not all mixed states 
like in the case of REE. Therefore it is usually easier to compute Eq rather than E R . We 
call the product state that achieves the maximum overlap with the closest product 
state (CPS). For general entangled states Eq gives the lower bound to E R but for pure 
stabiliser states the measures are equal [20] . 

Consider a pure state G %\ ® . . . <S>"Hjv of an -/V-partite system. It can be written 

as 

m = Z>i#>®--.®i^>, 

8=1 

where £j G C. The Schmidt measure [9] is then defined as 
Es(Wi) = logiCm, 

where R m %n is the minimal number of terms in the expansion of over all linear 
decompositions into product states. Note that Es is not continuous but this does not 
pose a problem as the set of graph states is also discrete. 

2.2. Graphs 

A graph is a pair G = (V, E) [21]. Elements of V — {1, . . . , iV} are the vertices and 
elements of E C [V] 2 are the edges connecting the vertices. Pictorially one represents 
graphs as a set of dots, representing the vertices, connected by lines according to E, 
representing the edges. We consider only simple graphs, that is graphs where the vertices 
are not connected to themselves so (a*, a*) ^ E, and the graph contains no multiple edges 
between the same set of vertices. Furthermore we concentrate only on connected graphs, 
that is any two vertices a^, aj G V are connected by a path in G. 

The neighbourhood of a vertex a G V , denoted by N a , is defined as the set of all 
vertices that are adjacent to vertex a, N a := {b G V\(a,b) G E}. The degree d(a) of a 
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Figure 1. Successive application of local complementation first to vertex 3 and then to 
vertex 1. The corresponding graph states of all three of these graphs are LC-cquivalent. 

vertex a is in our case the size of a's neighbourhood |iV a |. A vertex colouring is a map 
c : V —> S such that c(v) ^ c(w) when vertices v and w are adjacent. The elements of 
S are called the colours. A graph is called bipartite if it is two-colourable. If a graph 
cannot be coloured by only two colours then it is called non-bipartite. 

One crucial quantity that we make use of extensively is the independent set. An 
independent set is a set of vertices where no pair is adjacent. A maximum independent set 
is the largest independent set for a given graph and is denoted by a(G). A closely related 
concept is that of a vertex cover. A vertex cover is a set of vertices such that each edge 
of the graph is incident to at least one vertex in the vertex cover. Minimum vertex cover, 
(3(G), is then naturally the smallest such set of G. The relationship between a(G) and 
(3(G) is that they are complements of each other which means that a(G) + /3(G) = V. 
This means that finding one set automatically gives the other. However identifying 
either the maximum independent set or the minimum vertex cover is a known NP-hard 
problem [22] . 

Finally we mention a very useful transformation of a graph known as local 
complementation [23]. Local complement (LC) of G at vertex a, denoted by r a (G), is 
obtained by complementing the subgraph of G induced by the neighbourhood N a , and 
leaving the rest of the graph unchanged. It is equivalent to adding a fully connected 
graph of vertices N a , denoted G(N a ), to the original graph modulo 2 

r a :G^r a (G) :=G + G(N a ). 

Applying local complementation to a vertex a adds edges between its neighbours where 
the addition is performed modulo 2 as demonstrated in Fig. [TJ 



2.3. Graph states 

Consider a graph G and its associated graph state \G) which can be prepared by placing 
a qubit at each vertex in the state |+) = 4|(|0) + |1)) and applying the entangling 
control-Z gate CZy between all vertices % and j that are adjacent f2M I2"5] 

ig>= n cz t] \+r v . (i) 

An alternative and more efficient way of describing the graph states is using the 
stabiliser formalism. The graph state is the unique, common eigenvector in (C 2 ) v to 
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the set of independent commuting observables {gi}f =1 given by 

9i := X % (g) Zj. (2) 

The eigenvalues of the correlation operators in Eq. (j2j) are +1 for all i E V. The Abelian 
subgroup S of the Pauli group V v generated by the set of all the correlation operators 
{g,j\i G V} is called the stabiliser of \G). We can generate a basis for (C 2 ) y using the 
graph state \G) by defining the set of states 

9i (' A- ,.../.-,... /.\\ > = (— l) fc *|G f fc 1 ...fc i ...fcj V )> 
where k± . . . is a binary string. The original graph state given by Eq. (JTJ) is simply 
|Go...o)- These states are all locally equivalent as \Gk x ...k N ) = -^l^o.-.o) and 

therefore they all have the same amount of entanglement. The projector onto the graph 
state can be written in the following nice form 

l G >< G l = ^v5> ( 3 ) 

where the sum is taken over all elements of the stabiliser S. 

A graph state \G) corresponds uniquely to some graph G. However a situation 
may arise when two distinct graphs G and G' represent two graph states that are 
related by some unitary operation \G') = U\G). Their stabilisers transform accordingly 
S' = USU^ = {U(tW\o~ G S}. So it can easily happen that two seemingly different 
graphs represent graph states with the same amount of entanglement. In this paper 
we consider such unitary operations that permute the elements of the Pauli group and 
therefore map stabilisers to other stabilisers. These transformations are called local 
Clifford operations, C\ = {U G XJ(2) V \UP V W = V v }. So two graph states \G) and 
\G') are LC-equivalent iff they are related by some local Clifford unitary U G C\ . In 
fact the action of local Clifford unitaries on graph states can be described graphically 
using local complementation transformations on the corresponding graph [27] . Applying 
the local complementation transformation to vertex a of graph G yields a new graph 
r a (G). The corresponding graph states \G) and \r a (G)) are then LC-equivalent and are 
related by local Clifford operation \r a (G)) = U^\G) given by 



Two graph states \G) and \G') are LC-equivalent if their corresponding graphs are 
related by a sequence of local complementations G' = r 01 o . . . o r an {G). Furthermore 
the local Clifford unitary relating two equivalent graph states can be found efficiently 
[28]. 



2.4- Bounding the entanglement 



Finding a candidate separable state that can be used in quantifying the entanglement is 
often a very hard task. However proving that this state achieves the minimum value of 



Direct evaluation of pure graph state entanglement 



7 



relative entropy or the maximum overlap is equally difficult. Fortunately the situation 
for many classes of graph states is such that lower and upper bounds on the entanglement 
measures coincide therefore the exact value of entanglement is known. This makes direct 
evaluation of entanglement measures easier because it is enough to find separable state 
uj and product state |0) that yield this value for Er and Eq respectively. In the case of 
the Schmidt measure Eg, knowing its value tells us the minimum number of terms the 
linear decomposition of \G) into product states needs to contain. 

In [12] the authors showed how the upper and lower bounds for relative entropy 
of entanglement and geometric measure can be found. The upper bound is obtained 
by considering perfect LOCC discrimination of a subset of the complete orthogonal 
graph state basis {(G^...^)}. A lower bound on the number of these states that can be 
discriminated is given by maximising the number of stabiliser generators {g{\ that can be 
determined in a single setting of LOCC measurements. This bound can be achieved by 
identifying the maximum independent set a(G) of the graph G and measuring X on these 
qubits and Z on their neighbours. This can in fact be used to find an upper bound on 
the relative entropy of entanglement and the geometric measure N — \a(G)\ > Er = Eq 

Lower bound for the entanglement can be obtained in the usual way of relaxing 
the condition of full separability and maximising the entanglement E bi between all the 
bipartitions of the graph. This can be viewed as creating the maximum number of Bell 
pairs m p between the bipartitions by CZ and local Clifford operations applied within 
the bipartitions. Summarizing the lower and upper bounds for the relative entropy of 
entanglement and the geometric measure are [12] 

N - \a{G)\ >E R = E G > E bi = m p . 

Classes of graph states for which these bounds are equal include d- dimensional cluster 
states, GHZ states and ring states with even number of qubits. 

The bounds for Schmidt measure can be found in similar fashion [TBI EE]- The lower 
bound is given by the maximal Schmidt rank SR max which for graph states is equal to 
the entropy of entanglement SR max = —Tr[p A \ogp A ] where A C V is a subset of the 
vertices. The upper bound is given by the minimal number of local Pauli measurements 
needed to completely disentangle the graph state. This quantity is referred to as the 
Pauli persistency, PP(G), and is bounded from above by the size of the minimum vertex 
cover |/3 (Cr) | [21] • So the Schmidt measure is bounded as follows 

|/3(G)|>PP(G)>£ 5 >SR max . 

States up to 7 and 8 qubits for which these bounds are equal are categorised in [2E] and 
[T7] respectively. 

3. Evaluating entanglement measures 

Now we have all the necessary tools needed to construct the minimal linear 
decomposition of the graph state into product states and hence the closest separable 
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Figure 2. Star graph and complete graph are related by a single local 
complementation operation therefore their corresponding graph states are LC- 
cquivalent and in this case have entanglement of 1. Care needs to be taken when 
computing the lower and upper bounds as the entire LC orbit of a graph, not just the 
graph itself, needs to be considered. Red edges depict the maximum matching. 



state to and the closest product state But before we do this it is necessary to mention 
that in many cases the lower and upper bounds coincide. In this section we present an 
algorithm for constructing the closest separable states that saturate the upper bound. 
Therefore when the bounds coincide it is trivial to show that a candidate separable state 
is truly the closest one. However if the bounds are not equal then it is not possible to 
conclusively state that the obtained separable state is the closest one. It is therefore 
crucial to understand when the bounds are equal and when they are not. 

3.1. Equality of lower and upper bounds 

To understand when the lower and upper bounds are equal we first demonstrate how 
the lower bound for all three entanglement measures can be found in terms of a single 
quantity from graph theory. 

A matching M of a graph G is a set of independent edges. Two edges are 
independent when they do not have any common vertices. A subset of vertices U C V 
is called matched if every vertex a G U is incident to an edge in M. If the entire vertex 
set V is matched by M then it is called a perfect matching. 

The notion of matching is very useful because the cardinality of maximum matching 
l^fmaxl is equal to the maximum number of Bell pairs that can be created for a bipartition 
of G as well as the maximal Schmidt rank SR max . 

At this point we would like to stress one crucial point that has been implicitly 
covered in Section [2731 We have explained that two graph states \G) and \G') which are 
related by a local Clifford transformation have underlying graphs G and G 1 related by 
a series of local complementation operations. However the properties of the two graphs 
may be very different. In particular their maximum matching and minimum vertex cover 
may have different cardinalities. One such example is the GHZ state pictured in Fig. [2] 
This graph state corresponds to a star graph which can be transformed to a complete 
graph using local complementation. So the graph states corresponding to a star graph 
and a complete graph are LC-equivalent and therefore have the same entanglement of 
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1. However for the star graph we have |M max | = 1 and = 1 whereas the complete 

graph has \M max \ = [f J and |/9(G)| — N — 1. 

To overcome this apparent contradiction it is necessary to consider the entire LC- 
equivalency class of a graph G and compute the size of M max and f3(G) for all graphs in 
the equivalency class to find the minimum. Such canonical forms for all LC-equivalency 
classes have been identified up to 7 qubits in [161 126], up to 8 qubits in [17] and up to 
12 qubits in [29]. For a graph state up to 12 qubits it is enough to find which canonical 
form it is equivalent to. However for a general graph state of N qubits one has to 
first generate the entire LC orbit of the underlying graph. Therefore when computing 
the lower and upper bounds for the entanglement of a graph state it is not enough to 
consider only the graph itself. The entire LC orbit needs to be taken into account. In 
the rest of this paper we will assume that the canonical form of the underlying graph 
ith smallest M max and (3(G) across the LC orbit is known. 

Consider a graph G with a maximum matching M max . Select either endvertex 
from each independent edge in M max to form one partition A. The other partition 
is then given by all the other vertices in V, that is B = V — A. Apply control-Z 
gates locally within the partitions to remove any edges that do not cross the bipartition 
boundary. As a result of this only edges that contribute towards entanglement across 
the bipartition are kept. In certain cases this transformation leaves all the vertices in 
the partitions with either degree or degree 1. This corresponds to the case of Bell 
pairs being created across the bipartition. However in the general case some vertices 
will have degree higher than 1. Then we need to apply a series of local complementation 
transformations along with further control-Z gates to either decrease their degree to 1 
or completely disconnect them from the rest of the graph. Crucial observation here is 
that none of the above transformations delete an edge from the maximum matching 
M max . Therefore this procedure produces \M max \ Bell pairs across the bipartition. This 
method of transforming the graph state has been introduced in [12] and demonstrated 
for a selection of graph states with underlying bipartite graphs. Fig. [3] displays the 
procedure for a particular example of a non-bipartite graph. 

Having established that the lower bound for all three measures is given by the 
size of the maximum matching |M maa; | and the upper bound by the cardinality of the 
minimum vertex cover we will now investigate when these bounds are equal and 

when not. To do this we will partition the graph states into the set of bipartite and 
non-bipartite states, and consider them separately. 

Bipartite graph states. For all bipartite graph states the lower and upper bounds 
coincide. This can be seen immediately using Konig's Theorem (21] which states that 
for bipartite graphs the sizes of maximum matching and of minimum vertex cover are 
equal, \M max \ = \(3(G)\. 

Non-bipartite graph states. The situation for non-bipartite graph states is more 
complicated as Konig's Theorem does not hold anymore and in general |M maa; | ^ \/3(G)\. 
However there are still numerous cases when the two bounds are equal. In fact the crucial 
quantity that determines whether the bounds are equal is the size of the maximum 
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Figure 3. Non-bipartite graph (in [26] denoted as No. 19). The maximum matching 
M max = {(1, 2), (3, 4), (5, 6)} is represented by red colour. The bipartition is given 
by the sets A = {1,3,5} and B = {2,4,6}. Successive application of local control-Z 
gates within the partitions and local complementation results in a disconnected graph 
representing 3 Bell pairs. 

independent set |a(G)|: 

(i) If \a{G)\ < f then \M max \ ^ \/3(G)\. 

(ii) If \a(G)\ > f then \M max \ = \(3(G)\. 

(iii) If |a((j)| = y then we have the following two scenarios: 

(a) If M max is perfect then \M max \ = \(3(G)\. 

(b) If M max is not perfect then \M max \ ^ \(3{G)\. 

Statement 1 can be proven easily. For a graph state of N qubits the maximum 
number of Bell pairs that can be created is |M maa; | < [yj . However the upper bound 
for entanglement is given by N - \a(G)\ > f from our assumption. Therefore the 
bounds are not equal, |M maa; | ^ \/3(G)\. 

Statement 2 becomes clear when one considers a bipartition of the graph G with 
partition A given by the vertices in the maximum independent set A := {v\v G «(£?)} 
and partition B given by all the other vertices, that is the minimum vertex cover 
B := {v\v G (3(G)}. It is clear that all vertices v G A are not connected. Now consider 
the reduced neighbourhood of a vertex a G B defined as N^ ed := {v G A\(a, v) G E}. 
So is just the usual neighbourhood with vertices in B removed. When looking for 
the maximum matching of graph G we can consider each N^ ed separately. The crucial 
point to note here is that since all v G N^ ed are incident to a single vertex a G B 
we can only pick one edge that contributes to the maximum matching. The fact that 
l^marrl = follows straightaway. Including edges between vertices in B in the 

matching can only decrease the matching's cardinality so they have been omitted from 
the argument. This is illustrated in Fig. HI 

Statement 3 becomes clear after one realises what it means to have a perfect 



Direct evaluation of pure graph state entanglement 



11 



2 1 




by orange vetices. The reduced neighborhoods are N£ ed = {4}, N$ ed = {3} and 
j\fred _ {1^2}. The maximum matching is depicted by edges with solid lines across 
the bipartition A — B. 



matching. The upper bound for entanglement in both cases (a) and (b) is given by 
= N — \a(G)\ = y. If the matching M max is perfect then each vertex in the 
graph is matched. Since a matching is a set of independent edges this implies that 
l^maxl = y- Therefore the entanglement bounds coincide, |M maa ,| = On the 

other hand if M max is not a perfect matching then there exists at least one vertex which 
is not matched by M max . Hence the number of vertices that can are matched by M max is 
upper bounded by N — 1. We know that N is even so the size of the maximum matching 
is upper bounded by |M mcu ,| < ^y^. This shows that if |a(G)| = y and the maximum 
matching is not perfect then the bounds for entanglement do not coincide. 



3.2. Entanglement evaluation 

As mentioned in Section 12.31 the stabiliser S of an iV-qubit graph state generated 
by iV generators gi given by Eq. ([2]) stabilises a unique entangled state. Consider 
a scenario where we discard a set of k generators. The new stabiliser is given by 
iS>Ar_fc = (<7i, . . . ,gN-k) where we have relabeled the remaining generators from S for 
convenience. Because we no longer have a full set of N generators, Sw-k does not 
stabilise a unique state. Rather it stabilises a set of states {|^i), • • • , |^d)} that span a 
.D-dimensional subspace. The dimensionality of this subspace depends on the particular 
structure of the generators gi £ Sjv-fe and the states (l^;)} may be entangled or product 
states. We are interested in the case when <Sjv_fc stabilises a set of product states. 

Observation 1. (Stabilised entangled states): Su-h stabilises entangled states if it 
contains at least two generators g a ,gb £ <SN-k where g\, contains Z a . 

Such a scenario happens when the generators g a and gb correspond to adjacent 
qubits, that is (a, b) £ E. Generator g a stabilises some subspace % a spanned by states 
{lV , | a ' ) )} an d gb stabilises a subspace Tib spanned by states {IV'] 6 '')}- Since qubits a and 
b are adjacent, the action of g a on the states {|^ )} is to permute them and the same 
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is true for acting on )}■ In order for g a to stabilise the states )} we have 
to take superpositions of these states according to how they are permuted by g a . This 
finally yields states that are stabilised by S^-k, however due to the superpositions the 
stabilised states are entangled. We illustrate this in the following example. 

Example 1. Consider a three-qubit open linear graph state given by the stabiliser 
S = (XZI,ZXZ,IZX). 

91 92 93 

Say we discard the third generator so the new stabiliser is given by S 2 = {XZI, ZXZ). 
Each generator stabilises a set of states 

gi = XZI:{\+0.),\-l.)}, 

g 2 = ZXZ : {|0 + 0), |1 + 1), |0 - 1), |1 - 0}}, 

where the states of the form | + 0.) mean that the generator does not fix the last qubit. 
We can quickly check that the action of g\ on the states stabilised by g 2 is to permute 
them, that is <7i|0 + 0) = |1 — 0) and gi|l + 1) = |0 — 1). Requiring that g\ stabilises all 
the above states we are forced to take superpositions of the states which finally yields 
the states that are stabilised by S 2 to be 

{^(l°+> + I 1 ")) ® 1°)' ^d°-> + ® I 1 )}' 

The stabilised subspace is two-dimensional and because of the form of the generators 
of 1S2 it is spanned by entangled states. On the other hand we can quickly check that 
discarding g 2 from the original S produces a stabilised subspace spanned by product 
states {| +0+), I - 1-)}. 

The above observation tells us exactly when <Szv-fc stabilises a set of product states. 
The desired structure of the generators ^ G S^-h is achieved when the generators 
correspond to non- adjacent qubits. As a vertex colouring of the underlying graph G 
partitions the set of qubits into subsets of independent qubits, the possible ways of 
discarding generators from the original stabiliser S are given by all the possible vertex 
colourings that G admits. More specifically if the generators corresponding to qubits 
of the same colour are the only generators not discarded then SN-k stabilises a set of 
product states. Note that the converse is not true as it is possible for SN-k to contain 
generators corresponding to different colours and still stabilise a set of product states. 
This is illustrated in Fig. |5j 

Knowing how to discard generators from S to produce a product basis we can focus 
on the important question of doing this optimally. Ideally we would like to discard as 
few generators as possible. In light of the above discussion it is straightforward to 
see that the generators that need to be kept correspond to qubits in the maximum 
independent set a{G), or equivalently, the generators that need discarding correspond 
to the minimum vertex cover (3(G). Similar approach has been used to compute upper 
bounds for the three entanglement measures in [T21 [T6], [26] though our logic of deriving 
this result is complementary to the approach in these references. 
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Figure 5. Two valid colourings for a 5-qubit GHZ state, (a) Both 52 = (52,53) and 
$2 = (54, 55) stabilise a set of product states. Keeping certain generators corresponding 
to qubits with different colours can also stabilise a set of product states. Such a 
stabiliser is 1S3 = (.92,53,54)- (b) Maximum independent set, a(G) = {2,3,4,5}, 
identifying the largest set of generators that stabilise a set of product states. 



We will now show how the entanglement measures can be evaluated directly. We 
refer to the stabiliser of generators corresponding to qubits in the maximum independent 
set a(G) as S a , the stabilised subspace as % a and product states spanning this subspace 
as m)}. 

Theorem 1. (Minimal linear decomposition into product states): Given a graph state 
\G), its minimal linear decomposition into product states is a superposition of states 



where fi(S) is a binary- valued function and its value depends on the action of the 
original stabiliser S on the states l^f). This ensures that the form of \G) in Eq. @ is 
stabilised by the whole S. D a is the dimension of the subspace l-i a and depends on the 
size of the minimum vertex cover as 

Proof. The proof consists of two parts. First part shows that the decomposition in 
Eq. (jlj) actually gives the graph state \G). The second part shows that the decomposition 
achieves the minimum bound of the Schmidt measure Es(\G)). 

From Observation ([T]) we can be sure that the states \^f) are product states. 
Furthermore since these states are stabilised by S a we have g^f) = for all j G S a 
and i G {1, . . . , D a }. All we have to determine is the action of generators in the minimum 
vertex cover (3(G) on the states \ipf). Consider the action of one of the generators in 
the minimum vertex set g%, k G (3(G). Due to the structure of the correlation operators 
this permutes the elements in the stabilised set gk\ipf) = )■ If the qubit k is adjacent 
to another qubit in the minimum vertex cover then the action of may be to also 
introduce a negative sign, g^ipf) = — )• By taking superpositions of states |^f) that 
are permuted with each other and including the negative amplitudes we can construct 
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a new set of states |^-) that are stabilised by the generators ^ G a(G) as well as the 
new generator g^ G (3(G). The new states will be of the form 

Repeating this process with a new generator g t where I G (3(G) acting on the states |^) 
we obtain a new set of states stabilised by S a as well as g k ,gi- Extending this procedure 
to all of the generators in (3(G) we finally arrive at the state in Eq. (jlj). 

Showing that this decomposition is also minimal becomes now trivial for graph 
states for which |M maa; | = as log-Da reaches the lower bound for the Schmidt 

measure known from [16], [26]. Therefore the entanglement as evaluated by the Schmidt 
measure is given by 

E S {\G))= log D a , 

which concludes he proof. □ 

Observation 2. Bipartite graph states for which the maximum independent set a(G) 
is also given by a 2-colouring of G can be written as \G) = J2f=i IVC )• Because 
any qubit in the minimum vertex cover (3(G) has all its neighbours in the maximum 
independent set a(G), the action of any generator g k G Sp is to only permute the states 
stabilised by S a , that is g^ipf) = IV'f) f° r hi e a(G). Therefore all the amplitudes in 
the linear decomposition of Eq. @ are positive. 

Using the decomposition in Eq. (jlj) we can now compute the other two measures. 

Theorem 2. (Relative entropy of entanglement): The closest separable state u to a 
given graph state \G) is given by 

Ua i=i 

So the closest separable state is given by equal mixture of the states stabilised by S a . 

Proof. States of the form given by Eq. (JSJ) are clearly separable because they are a 
mixture of product states. Computing the relative entropy between p = \G)(G\ and u 
we get 

S(p\\u) = - Tr[plogu;] 

ijk 

= logD a . 

For graph states which have |M maa; | = this is the lower bound found in [T2] . 

Therefore u is the closest separable state and the relative entropy of entanglement is 
given by 

E R (\G)) = log D a . 

This concludes the proof. □ 
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The closest separable state u can be expressed in a similar form to Eq. (j3J) as the 
sum of all the elements in S a . 

Observation 3. We can write the closest separable state as 

This can be most easily seen by using Eq. §5§ and noting that (^- Y^i=i l^f) (^f I) iV'St) 
Using Eq. (E} and the fact that cr^f ) = |^f) for all a G «S Q we can obtain 

2l«(G)l 



1 __ 2\ a(U " 



2Ij9(G)I 

where we used |a(G)| + = N and that the cardinality of <S a is 2' a ^ G ^. This 

establishes equivalence between the two forms of u in Eq. (JSJ) and Eq. fl6]). 

Theorem 3. (Geometric measure): The closest product state \4>) is given by any state 
from the stabilised set {l^f )}, 

|0) = IC), Vze{l,...,D a }. (7) 

Proof. Using the minimum linear decomposition into product states in Eq. (j3J) and 
substituting \<f>) = \ipf) we find that the geometric measure is 

E G (\G})= -l g|^f|^=X;(-l)^^|^)| 2 



3=1 

= log A*- 

This is the lower bound found in [T2] and therefore concludes the proof. □ 



We demonstrate how to quantify the entanglement measures in the following 
example of a 6-qubit graph state. 

Example 2. Consider a graph state \G) of 6 qubits with an underlying graph G pictured 
in Fig. [6j The graph state \G) is stabilised by 

S = ( XIIIIZ , 1 XIII Z , IIXIZI , 

9i 92 93 

iiixzi , IIZZXZ , ZZIIZX) . 

94, 95 96 

The maximum independent set can be quickly obtained by a 3-co louring and is a(G) = 
{1,2,3,4}. Therefore the corresponding stabiliser is S a = (g±, g 2 , g$, #4). |«(G)| > 3 
which means that the lower and upper bound coincide. The stabilised states given by 
S a are 

{!#)} = {! + + + +00), I - - + +01), 1 + + - -10), 1 11)}. 
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Figure 6. Graph G considered in Ex. ([2]). (a) Graph G is bipartite since a 2-colouring 
can be found, (b) However the maximum independent set a(G) is identified using a 
3-colouring. 

The action of the stabiliser Sp associated with the minimum vertex cover is the following 

= m and g 8 \W) = -|^), 
g 6 \^) = m) and g e \^) = -|<>, 
which means that the graph state can be written in the following form 

\G) = \[m+m + m-m)\. 

Using the closest separable state u in Eq. (JSJ) and the closest product state in Eq. flTj) 
we can finally show that for this graph state \G) 

E S (\G))=E R (\G)) = E G (\G))=2. 

Knowing how to find the closest separable state to to a given graph state \G), it is 
now possible to determine the form of the closest separable states to all of the graph 
states in the orbit generated by local Clifford operations. Two graph states are LC- 
equivalent, \G') = U LC \G), if they are related by a local unitary U LC . Using Eq. (]1J), 
\G') can be expressed as a superposition of states U LC \il)f). Applying the unitary we 
see that the closest separable state has the form u>' — ^ U (ipf l \U LC ^ . 

3. 3. Maximal entanglement of graph states 

An immediate consequence of this approach to analysing entanglement is that we can 
identify which graph states are maximally entangled in many cases. By maximally 
entangled we mean with respect to the three measures we are considering. For any 
graph state that is maximally entangled it must be true that Er = Eq > \M max \ and 
Es > |-Wmax|- F° r bipartite graph states these measures cannot be larger than this value 
due to Konig's Theorem [21]. Therefore any bipartite graph state \G) whose underlying 
graph G has the property that \M max \ = [yj is maximally entangled. Examples of such 
states include linear graph states with open boundaries, ring states with even iV and 
cluster states in d dimensions. This holds true also for non-bipartite graph states for 
which the bounds are equal. An example of such a state is pictured in Fig. HI 
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Figure 7. Four lattices for which we identify the maximum independent set a(G) 
allowing us to determine the scaling of the gap between upper and lower bound for 
entanglement for general N. The maximum independent set corresponds to orange 
vertices. The gap is finite for the triangular lattice (a), the kagome lattice (b) and the 
hexa-triangular lattice (c). For the hexagonal lattice (d) it is trivially zero since this 
lattice is bipartite. 



The situation is quite different for the case of graph states with unequal 
entanglement bounds. Again any graph state that is maximally entangled must have 
l^maxl = LtJ • However the true value of entanglement is now unclear. For general 
graph states it is not even possible to calculate this bound. Therefore we limit ourselves 
to various types of regular two dimensional lattices for the rest of this section. 

The lattices that we consider are pictured in Fig. [7J They are the triangular, 
kagome, hexa-triangular and hexagonal lattices. The triangular, kagome and hexagonal 
lattices have been shown to be universal resources for measurement-based quantum 
computation [30]. We are interested in the scaling of the gap between upper and lower 



bounds A = — |M max |. For hexagonal lattice this gap is trivial since the lattice 



is bipartite and therefore A 



hexagonal 



Hri 



*hex—tri 



N' 

y 



0. For the other three lattices that gap is 
E (^-3j), 



N-2 



(12JV - 3 V9 + 12N + 9) - 
(6iV- V13 + 3JV-11) 



N 

y 



N 

y 



where A tri is valid for L > 3, with L 2 = N being the number of vertices on one side of 
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Figure 8. Scaling of the gap A for the three lattices pictured in Fig. [Bfa)-(c). The 
gap A increases linearly in the leading term. 



the triangular lattice. The particular form of A depends on the boundary conditions. 
However the general behavior of the scaling remains unchanged for different boundary 
conditions. Interestingly the lower bound for all four lattices is the same, |M maa .| = [yj . 
The scaling of A is pictured in Fig. |HJ 



4. Alternative description of CSS 

In this section we focus only on one measure, the relative entropy of entanglement, and 
we abandon the stabiliser generator description of graph states to see if it is possible 
to arrive at the closest separable state u using complementary methods. We consider 
two approaches. The first one is inspired by projected entangled pairs states description 
of graph states. Usually PEPs methods are used to describe pure entangled states, 
however we adapt this approach to construct the closest separable state u. The second 
approach of obtaining u relies on introducing optimal amounts of noise in the form of 
relative phases and averaging over these phases. The success of both of these methods 
rests on our ability to identify the maximum independent set a(G). 



4-.1. PEPs construction 

We briefly highlight the PEPs description of entangled states [IB] . Consider a graph 
state \G) with iV qubits. A PEPs |^) G C 2 " is constructed by replacing a physical qubit 
a by |iV | virtual qubits where |7V | denotes the degree of vertex a. Each physical edge 
(a, b) is then replaced by maximally entangled state of the corresponding two virtual 
qubits \G 2 ) = 4s(|0+) + |1— )). The original graph state \G) can then be obtained by 
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Figure 9. (a) The two 2-qubit states from Eq. |8] used in construction of closest 
separable states. The dashed edges represent that the states are separable. We choose 
the convention that orange vertex corresponds to a qubit in {|±)} basis and blue vertex 
corresponds to a qubit in {|0), |1)} basis, (b) Open linear 4-qubit graph state IG4) and 
its corresponding closest separable state W4. 



applying a projector P a : = |0)(0| al . . . (0| alJVal + |l)(l| al . . . (l| a|JVa at each physical site. 

The above approach can be adapted to describe separable mixed states with few 

changes. Instead of maximally entangled pairs of virtual qubits the basic building blocks 

are maximally correlated separable pairs of qubits u?. Maximally correlated in this sense 

means that the relative entropy between the separable states and the tensor product of 

its subsystems is unity, S^o^Ua^ ) = S(uj2\\jl ® I) = 1 where u4 = Tr 2 [^2] and 

(2) 

similarly for uj 2 . In fact it will be necessary to use two various separable states uj^. The 
projectors being applied to physical sites will also have a different structure compared 
to the case of pure entangled states. 

Define two 2-qubit maximally correlated separable states of virtual qubits i' and f 

<4 r :=|+0)(+0| + |-l)(-l| 

^:=|0+)(0+| + |l-)(l-|. 

Here and in the rest of this subsection we omit the normalisation constants. Virtual 
qubits are denoted by primed letters a' and physical sites by a. These states are in 
fact both closest separable states to a 2-qubit graph state. It is also useful to give these 
states a graphical representation depicted in Fig. HI Using the separable states in Eq. (jSJ) 
we construct a separable state of virtual qubits by placing either ui^-, or ui^-, on edges of 
the graph state. The 2-qubit separable states are picked in such a way that all virtual 
qubits at a physical site a G a(G) are orange. In the case of bipartite graph states this 
means that all virtual qubits at physical sites b G (3(G) will be of the same colour, blue. 
However this is not true anymore in the case of non-bipartite graph states where virtual 
qubits at a physical site b G (3{G) will be of both colours. This is illustrated in Fig. EJ 
Finally in analogy to the usual PEPs construction the virtual qubits at physical sites 
are projected using the following maps 

P a A := |0)(0...0| + |1)(1...1|, 

H+X+l + l-X-l- 
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The projector is applied if the physical site a G (3(G) contains at least one virtual 
qubit of blue color and P B is applied if a G a(G) which means that all its virtual qubits 
are orange. |+) is an equal superposition of all tensor product states {|+), |— )} which 
are +1 eigenstates of X ® . . . <g) X. Similarly |— ) is an equal superposition of all — 1 
eigenstates of X <S> ■ ■ ■ (8> X. 

How this construction works is most easily seen when illustrated on an explicit 
example. Consider an open linear graph state of 4 qubits as in Fig. [HJ Following the 
2-colouring of the graph the tensor product of 6 virtual qubits takes the following form: 

The desired 4-qubit closest separable state is then obtained by applying the projector 
u, = (P 2 A ®P 3 B )cu' 6 (P 2 A ®P 3 B ), 

where P 2 A = |0)(00| + |1)(11| acts on physical site 2 and P 3 B = ]+)(+ + | + |+)( | + 

l — ) ("I — I + l — )( — I" I ac ^ s on physical site 3. 



4-2. Noise construction 

In this subsection we focus on the third approach of constructing the closest separable 
state u. The common feature of both previous approaches is the necessity of identifying 
either the maximum independent set a(G) or the minimum vertex cover (3(G). This 
remains true for the approach presented here as well. However this time we ask the 
question if there is a simple way of obtaining u by introducing noise to the pure state 
rather than resorting to methods based on stabiliser generators or pairs of maximally 
correlated separable states. 

This new approach relies on introducing distinct relative phases to certain qubits 
and then averaging over them to obtain uj. The minimum number of the relative phases 
is equal to the cardinality of the minimum vertex cover 

The graph state vector can be written in the following form [24] 

1 N 

|G> = 2i^®(|0>i + |l>i^v}). (9) 

3=1 

where the Pauli Z matrix is applied to a subset of qubit j's neighborhood TVj := {i G 
Nj\i > j}. Now lets define a new vector |$) that differs from the graph state in Eq. (J9]) 
in that it contains the above mentioned relative phases <pj 

1 N 

m--=7^- 2 ^)(\0h + e^ mU) \lhZ Nf X 



2 N/2 

i=i 

where m(j) : V — > {0, 1} is an indicator function from the set of all vertices V given by 

if jea(G), 

1 if j G /3(G). 

So the relative phase <pj is applied only to qubits that correspond to vertices in the 
minimum vertex cover. For example a 3-qubit open linear graph state with a relative 



m(j) :-- 



Direct evaluation of pure graph state entanglement 



21 



phase is |$) = ^ 5 (|0) 1 + |l) 1 Z 2 )®(|0) 2 + e^|l) 2 Z 3 )(2)(|0)3 + |l>3). Finally the closest 
separable state is given by averaging over these phases 

where d<p = n {i | je(3(G)} c%. 
5. Conclusions 

We have presented a method of evaluating three multipartite entanglement measures in 
pure graph states. Our approach uses simple group theoretic arguments to identify a 
suitable subspace of the original Hilbert space whose properties can be used to find the 
relevant closest separable and product states as well as the minimal linear decomposition 
of a pure graph state |G). 

The problem of evaluating entanglement measures can be mapped directly to a well 
known problem of identifying the maximum independent set in graph theory. Knowing 
the size of the maximum independent set corresponds to knowing the upper bound for 
all three entanglement measures as well as whether the upper and lower bounds are 
equal. Identifying which qubits comprise the maximum independent set allows us to 
construct the minimal linear decomposition of \G) into product states as well as its 
closest separable and closest product state. 

The closest separable state u admits a non-stabiliser description using a PEPs 
inspired construction. This immediately begs the question whether a suitable u can be 
constructed for weighted graph states. Any realistic scheme of preparing graph states 
will use entangling gates best described by a general control phase gate where the phase 
does not always have the ideal value of <fi = ir. Rather it is picked randomly from some 
distribution centered around this ideal value. 

Our techniques developed so far are the first step towards investigating total 
entanglement properties of such realistic systems. Furthermore our methods can be 
used in the study of entanglement in lattice models with long-range interactions. The 
strength and range of these Ising-type interactions can be captured by the phase in the 
entangling gate. All work on this topic has so far been limited to the study of bipartite 
entanglement measures. We will present our findings on this topic in a separate paper. 

Interesting observation is that the lower bound for entanglement of any pure graph 
state can be found efficiently since the maximum matching problem of an arbitrary 
graph is in the P complexity class. On the other hand the maximum independent 
set problem, and hence also evaluation of the upper bound, is NP-hard for a general 
graph. One exception are bipartite graphs for which |M max | = N — |a(G)|. Efficient 
estimation of the upper bound is still an open problem. It would be interesting to see 
whether separable states that approximate the closest separable state can be constructed 
efficiently. 

A closely related question is whether there is some deeper relationship between the 
upper and the lower bound. Calculations for graph states up to 10 qubits suggest that 
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the gap between the bounds grows very slowly. It would be interesting to see how this 
changes for the case of larger and more general graph states. We have demonstrated 
that for certain regular lattices in two spatial dimensions the gap A increases linearly 
with the number of qubits N. This is not too surprising. Because we have considered 
regular lattices, linear increase of N results in a linear increase of the cardinality of 
the maximum independent set |a(G)| while the size of the maximum matching |M max | 
remains constant. However this behaviour is not likely to be true for more general 
graphs. Particularly interesting would be to find a relationship between the size of 
the gap and some structural quantities of the underlying graph that can be computed 
efficiently 

This naturally leads to the final question which is concerned with evaluating 
entanglement in graph states where the bounds are not equal. In this case our methods 
can achieve the upper bound and therefore do not say anything concrete about the 
actual entanglement of the graph state. Numerical evidence suggests that for certain 
states geometric measure is less than the upper bound [31] whereas for some other states 
it is equal to the upper bound. An open question is to see if the three entanglement 
measures are still equal when the bounds are different and to characterize states whose 
upper bound is the actual value for entanglement. 
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